Continuous monitoring can improve single-photon probability 
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An engineering technique using continuous quantum measurement together with a change detec- 
tion algorithm is proposed to improve the probability of single photon emission for a quantum-dot 
based single-photon source. The technique involves continuous monitoring of the emitter, integrat- 
ing the measured signal, and a simple change detection circuit to decide when to stop pumping. The 
idea is to pump just long enough such that the emitter + cavity system is in a state that can emit 
at most one photon with high probability. Continuous monitoring provides partial information on 
the state of the emitter. This technique is useful when the system is operating in the weak coupling 
regime, and the rate of pumping is smaller than, or comparable to, the emitter-cavity coupling 
strength, as can be the case for electrical pumping. 
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I. INTRODUCTION 

Generation of single-photon states has wide-ranging 
applications, spanning quantum computing, quantum 
imaging, metrology, communication, and cryptography, 
amongst others [3]. These applications are important to 
scientific and technological progress in many important 
areas. For the current work, we are especially interested 
in Linear Optical Quantum Computation (LOQC) [Ul2]- 
A key requirement of LOQC is the availability of high 
quality single-photon states on demand. 

Semiconductor quantum dot-based implementations of 
single-photon sources are of particular interest, as they 
scale well upon integration and are amenable to commer- 
cial fabrication techniques [3JH]. Typically, implementa- 
tion of these devices involves a quantum dot (QD or dot) 
inside a microcavity, with the dot acting as the photon 
source: the presence of the cavity increases the collection 
efficiency due to spatial confinement of the photons [3H5] . 
A wide variety of cavities with different sizes, shapes, and 
quality factors have been designed and fabricated [5]. 

A quantum dot in a microcavity can be pumped ei- 
ther optically or electrically. Optical pumping is more 
straightforward experimentally, but electrical pumping 
may be better suited to large scale integration. Also, 
as it does not directly insert photons into the cavity, it 
opens up the possibility of pumping directly into an en- 
ergy level resonant with the cavity mode, which may re- 
duce timing uncertainties in photon emission [6]. It also 
allows a channel for continuous measurement of the dot. 
In electrical pumping, a bias voltage is applied across a 
quantum dot p — n diode that enables an electron to tun- 
nel through from n-type onto the dot (present at p — n 
junction). Once an electron tunnels through to the dot, 
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no new electron can tunnel through due to the Coulomb 
blockade effect [7 . A further increase in the bias voltage 
then enables a hole from the p-type to tunnel to the dot. 
Recombination of electron and hole in the dot follows, 
leading to a photon emission into the cavity. Recombi- 
nation also leads to a drop in the potential across the 
quantum dot p — n diode. By observing this change in 
the potential across the quantum dot p — n diode, we can 
gain information about the state of the dot. This obser- 
vation can be seen as a weak quantum measurement that 
provides partial information about the system [8]. 

A good single-photon source should be able to pro- 
duce exactly one photon at the required time, in a 
specified state. Indistinguishability is a measure that 
captures the specificity of the photon state (due, e.g., 
to time-uncertainty of photon emission), while single- 
photon probability determines the likelihood that a single 
photon is indeed emitted. We have considered continuous 
monitoring as a tool to improve indistinguishability else- 
where [6]. We focus on single-photon probability here. 

There are various processes that affect single-photon 
probability: collection efficiency (emitting into the cav- 
ity mode), pumping, cavity leakage to non- waveguide 
modes, and photon loss, to name a few. Many of these 
processes depend on the fabrication techniques and the 
materials that go with it, and hence can only be im- 
proved by building better sources. Some, however, may 
be improved by better control. The most obvious pro- 
cess to treat from this point of view is pumping: how 
long should we pump at a given pumping strength to 
maximize the single-photon probability, given that other 
parameters are fixed? 

In the case of strong pumping — as optical pumping of- 
ten is — the answer is straight forward: pump for a dura- 
tion short compared to the emission rate of the dot. How- 
ever, weakly pumped systems (like electrical pumping) 
present a more complicated scenario. Since the pumping 
strength is weak, we might have to pump for times com- 
parable to, or even longer than, the emission rate of the 
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dot. We thus may increase the likelihood of multi-photon 
emission. Since electrical pumping with semi-conductor 
QD based implementations is of active experimental in- 
terest [7 j, understanding the behavior of weakly pumped 
systems is critical. 

In this work, we use continuous quantum measure- 
ments to improve single-photon probability. The idea is 
simple: to monitor the state of the emitter continuously, 
and determine when to stop pumping energy into the sys- 
tem based on the information obtained. Unfortunately, 
the output from monitoring such a microscopic system 
is intrinsically noisy, and little time is available to pro- 
cess the received signal. This complicates the procedure. 
We utilize a sequential statistical technique called Cu- 
mulative SUM (CUSUM) as the decision making process. 
We show numerically that this mechanism substantially 
improves single-photon probability in the weak coupling 
regime, particularly when the pumping rate is compara- 
ble to the QD-Cavity coupling strength. 



A. Overview of the paper 

The emitter we consider is a p-n diode operated as 
a single-photon LED, though the techniques we discuss 
are probably applicable to other systems. We discuss the 
emitter in Sec. [TTJ We then capture the essential features 
of the LED's operation in a schematic model in Sec. |III| 
we present details regarding the system in a cavity-QED 
setting, and give its energy-level diagram; we present a 
stochastic master equation including the relevant phys- 
ical processes, and discuss the parameter regime of op- 



eration. Section [TV] presents the change detection algo- 
rithm and the decision circuit used to stop pumping; we 
consider sequential CUSUM technique in Sec. |IV A| and 
present Bayesian approach in Sec. IV B We give results in 
Sec.[V| we discuss the deterministic case first in Sec. |VB[ 
before proceeding to analyze the numerical performance 
of the CUSUM-based technique in Sec. |V C| Finally, we 
scrutinize the effect of monitoring efficiency on perfor- 
mance in Sec. IV Dl 



II. SINGLE-PHOTON LIGHT EMITTING 
DIODE 

Figure [I] shows a quantum dot p-n diode that acts as a 
single-photon source. A dot is present at the p-n junction 
and is assumed to be inside an optical microcavity in 
the weak coupling regime. We assume that the diode is 
forward biased and is in the Coulomb blockade regime. 
When biased at an appropriate level, an electron (e~) 
tunnels through to the dot from the n-side; this electron 
remains in the dot until a hole (e + ) tunnels through from 
the p-side. This leads to recombination of e~ and e + , 
and a photon is emitted into the cavity. The photon 
subsequently leaks out to an external mode, such as an 
outgoing waveguide. 




Quantum dot embedded in a cavity 



FIG. 1. (Color online) Quantum dot p-n diode or p-i-n 
hetero junction. A quantum dot p-n diode comprises an 
insulator sandwiched between p- and n-type silicon. A quan- 
tum dot is fabricated inside the insulator, and this is con- 
tained within an optical microcavity. The diode is biased in 
the forward direction, such that a single electron (e - ) tunnels 
through from the n-side to the dot. The electron remains in 
the dot until a hole (e + ) tunnels through to the dot from the 
£>-side. The electron- hole pair in the dot recombines to emit 
a photon into the cavity, which subsequently leaks out to an 
external mode. 



While photon generation as described above is intu- 
itively straightforward, we need to control the pump- 
ing so that the system generates at most one photon 
with high probability. Since the dot is inside a micro- 
cavity, there is a non-zero probability of multi-photon 
emission; this happens when there is more than one re- 
combination event in a pumping cycle. Determining the 
length of a pumping cycle should ideally be determined 
by the knowledge of when the first e~ tunneling event oc- 
curs. This is difficult to know, as the tunneling process is 
stochastic. In the absence of specific information about 
tunneling times, the best that can be done is to time 
the pumping a priori, either to maximize the probability 
of a single photon, or keep the multi-photon probability 
below a given threshold. (These are not necessarily the 
same thing.) In principle, though, we can do better. 

We continuously monitor the state of the dot by mea- 
suring the voltage across the hetero junction; the output 
record gives information about whether an e~ has tun- 
neled onto the dot. We use this information to stop 
the pumping cycle when the dot is in its excited state 
(equivalently, after an e~ tunneling event). If done cor- 
rectly, the dot is excited only once, and the diode emits 
at most one photon. However, the measurement signal 
is noisy, and the decision of whether the dot is in its ex- 
cited state or not is not necessarily easy. We use a well 
known sequential statistical decision technique known as 
Cumulative SUM (CUSUM), a decision circuit with a 
simple imple men tation that accommodates noisy obser- 
vations (Sec. IV). 
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III. SYSTEM MODELING AND PARAMETER 
REGIME 

Figure [2] schematically models the single-photon LED 
operation described in Sec. [TTJ The system comprises 
three quantum objects: a quantum dot, an optical mi- 
crocavity, and an external mode (e.g., a waveguide). We 
assume that the dot is electrically pumped: this allows an 
electron to tunnel through from the n-side to the dot, and 
recombine with a hole from the p-side, when the diode 
bias voltage is favorable. The bias voltage is controlled 
externally, giving us some control over the e~ tunneling 
event. To exert this control effectively, we continuously 
monitor the state of the dot and use the output to decide 
when to stop pumping. 

We include the external mode in our description 
in order to calculate the various photon emission 
probabilities — p(0), etc. These probabilities are 

the quantitative measure of "goodness" of the system as 
a single-photon source. Obviously, we want the single- 
photon emission probability, p(l), as close to unity as 
possible; unfortunately, in many parameter regimes p(l) 
is not very close to 1, especially if it is electrically 
pumped. In this case we must explore a more com- 
plex trade-off landscape, as it is no longer sufficient to 
just maximize p(l); we must also consider how high a 
zero-photon probability p(0) and multi-photon probabil- 
ity p{2 J r) we can tolerate. For certain applications it is 
tolerable for the single-photon source sometimes to emit 
no photon, but multi-photon emission must be strongly 
suppressed. Thus, in our trade-off analysis we impose an 
additional constraint on p(2+), and try to maximize p(l) 
subject to this constraint (see Sec. |V|). 

A. The system 

The dot in our model has 2 energy levels: the ground 
state |G), and the excited state \X). The cavity and 
external mode contain some number of photons, and are 
represented by the usual photon number state notation — 
|0), |1), 1 2) etc. Thus, any state in this system has 
the form \G/X, 0/1/2..., 0/1/2...) where the order cor- 
responds to the dot, the cavity and the external mode, 
respectively. 

We assume that the dot's excited state \X) is reso- 
nantly coupled to the cavity mode. The interaction be- 



T-L and V are superoperators: 

H[A]p = ApA ] - (A ] Ap + pi f i)/2, 



tween the dot and the cavity is given by the Jaynes- 
Cummings Hamiltonian: Hi = ihg{o)a~ — a<r + ), 
where g is the interaction strength, a) (a) is the cre- 
ation (annihilation) operator acting on the cavity mode, 
and a~ — |G)(X|; we operate in the interaction picture, 
and henceforth set the total Hamiltonian H to Hj. The 
system is initially decoupled, in the state |G, 0,0). 

Incoherent processes included in this model are pump- 
ing, spontaneous decay of the quantum dot, cavity leak- 
age and dephasing. Q is the rate of pumping; T is the 
spontaneous emission rate for the X G transition; k is 
the photon leakage rate from the cavity to the external 
mode, and 7 is the dephasing rate between emitter and 
the cavity mode. 

We treat electrical pumping as a incoherent pro- 
cess [12 . The process of pumping involves an electron 
tunneling through to the dot when the bias voltage at the 
diode junction is favorable. Since electron hopping (re- 
side —> dot) happens at random times, we model pump- 
ing as an incoherent G — >• X transition with rate Q. 

Figure |3] shows the energy level diagram of the single- 
photon LED system. The system moves up the en- 
ergy ladder as we pump longer. For finite pumping 
time, the system eventually evolves to its possible final 
states I G, 0, 0/1/2 ... ), with corresponding probabilities 
p(0/l/2 . . . ). Both the cavity mode and the external pho- 
ton mode in principle have infinitely-many energy levels; 
however, since our goal is to generate a single-photon 
state |G, 0,1) with high probability (p(l) — » 1), and we 
are assuming that the coupling between the cavity and 
external mode is stronger than the coupling between the 
dot and the cavity, we truncate the higher energy states 
(> 2) as shown in Fig. [3] By truncating the state space, 
we make the state |G, 0, 2) now represent |G, 0, 2+) — that 
is, a state with 2 or more photons in it — and its corre- 
sponding probability becomes the multi-photon probabil- 
ity p{2 J r). Note that this approximation neglects some 
effects that could in principle contribute to the single- 
photon probability, such as a series of reabsorptions and 
spontaneous emissions. Since the probability of the sys- 
tem being in the higher energy states is small in the first 
place, and the spontaneous emission rate is assumed to 
be quite low, the neglected effects should have little im- 
pact on the estimates of single-photon probability. 

The evolution of the system is described by a stochastic 
master equation (SME) [TP] : 



(i) 

1 

V[T X ]P = V x p + pV x - 2{P X )P, 
where (Vx) = Tr{Vxp} is the quantum expectation; 



dp =- % -[H,p]dt + (TH[&-] +nU[ab ] } + iH[P x ] + nn[a + }^jpdt+ y/rnT>[Px]pdWt. 
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FIG. 2. (Color online) Single-photon LED operation: Schematic model. A single-photon LED consists of three quantum 
objects: an emitter (QD), a microcavity and an external mode. Electrical pumping involves tunneling of an electron (e - ) from 
the n-side to the dot; an e~ — e + recombination leads to an emission of a photon into the cavity, which then leaks out to the 
external mode. The e~ tunneling occurs when the bias voltage in the diode is favorable to the event. We control the bias 
voltage externally; we keep the bias 'on' to allow e~ tunneling, and turn it 'off' to stop tunneling. We monitor the state of the 
dot continuously, and use a sequential algorithm called CUSUM and its associated decision circuit to decide whether to keep 
the bias 'on' or to turn it 'off.' 



dW t is a Brownian motion with 

E[dW t ] = and E[dW t dW s ] = 5(t-s)dsdt, 

which characterizes a Wiener process; here E is the ex- 
pectation of a random variable. Also, <j + = |X)(G|, is 
the creation operator acting on the external mode, and 
Vx = \X)(X\ is the dephasing or "observer" operator 
acting on the quantum dot. Dephasing has two potential 
sources: the interaction of the emitter with other de- 
grees of freedom, for instance phonon modes in the dot 
[5], and the back action of a measuring device coupled 
to the dot. A measuring device will allow us to recover 
some information about the system of interest, but it is 
unlikely that we can tap into internal modes of the dot; 
we use < rj < 1 to denote the measurement efficiency 
with which we (the "observer") recover information lost 
in dephasing. 



B. Parameter regime 

We operate in the weak coupling regime, and assume 
that k is the dominant system parameter. We require 
that spontaneous emission T be small, for higher T im- 
plies lower single-photon probability. 



The technique developed in this work is based on con- 
tinuous monitoring of the dot, and will be useful only 
when the pumping rate Q is comparable to g. This is 
because it takes time to gather information using contin- 
uous monitoring, and with higher Q the chances are that 
the decision to turn off the pumping will be too late, 
increasing the multi-photon probability above tolerable 
levels. In fact, the scenario where g « ft has a much 
simpler solution: turn on pumping for time much less 
than the emission time scale (~ (g 2 /^) -1 ); this will work 
quite well because we move up the energy ladder (Fig. [3| 
only when the dot makes multiple X — »• G transitions; 
since T is very small, the transition time scale is thus set 
by g- 

We therefore are interested in the case where the pa- 
rameters satisfy the following conditions: 

r<<^,#<7<^, (2) 

for our technique to be useful. 

IV. CHANGE DETECTION ALGORITHM 

The output signal obtained from our continuous mea- 
surement is given (in rescaled units) by 

J(t) = (V x )(t) + 0£(t) (3) 
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Higher energy states (3+ quantum) 




FIG. 3. (Color online) Energy level diagram. A state of the system is described by three quantum numbers: the dot, the 
cavity and the external mode. In this diagram, we show those states with at most 2 quanta of energy and their significant 
dynamical processes: broken arrows indicate incoherent processes, and solid arrows indicate coherent evolution. If the pumping 
duration is finite, \G, 0, 0), \G, 0, 1) and \G, 0, 2) are the possible final states of this system, and their corresponding probabilities 
represent the zero photon, 1 photon and 2 photon probabilities, respectively. 



where f3 = (^7) -1 ^ 2 and = dW t /dt is Gaussian 
white noise with zero mean, i.e. K[dWt] = and dW 2 = 
dt, where E is the expectation of a random variable 
[TO] : Vx = \X)(X\ is the dot's excited state projec- 
tion operator, while (Vx)(t) is its quantum expectation 
Tr{V x p(t)}. 

The excited state of the dot is |X), and the unexcited 
state is \G). The measurement output conveys the infor- 
mation about this state: 

mn-AAi /A f0 if the dot is in |G), , AS 

nj(t)] = {Vxm = 1 1 ifthedotisin |4 (4) 

as E[£(t)] = 0. If the output contained no noise (/? — >• 0) 
then the state of the dot would tend to be localized at ei- 
ther \X) or |G), and transitions between would be readily 
detectable. The presence of noise complicates matters, 
since it can mask the state of the dot — indeed, if the 
noise is high, it is easy to miss the transitions. Note that 
if we average the output signal over a short interval At, 
its variance is a constant: a 2 = (3 2 /At. When an elec- 
tron tunnels through, there is a change in the mean of 
the output signal. By detecting this change, we could de- 
tect when a tunneling event takes place, and turn off the 
pumping. However, we cannot access the mean E[J(t)] 
directly, but have to infer it from J(t) [Eq. 

We obtain the output signal J(t) continuously in time, 
and the decision to turn off the pumping must be made in 



real time. This imposes practical restrictions on the kind 
of algorithms that are feasible; for instance, Bayesian 
machine learning-type algorithms [13] may be too slow 
to be useful in practice. Sequential algorithms [13] are 
procedures that use only the output signals gathered to 
the present time, and not a priori information, to infer 
the probability density function (pdf). As such, they are 
generally suboptimal, and will succeed only in certain pa- 
rameter regimes, but they have the virtue of being easier 
to implement in practice. One such algorithm is Cumu- 
lative Sum (CUSUM). 



A. Sequential CUSUM procedure 

We present a discrete version of CUSUM here for sim- 
plicity, which can straightforwardly be extended to the 
continuous case. Let y\ = J(£i),2/2 = Jfe),... repre- 
sent the time series of the output signal, averaged over 
a succession of intervals of size At. We know that Yk 
is a 2 parameter Gaussian random variable with a vari- 
able mean (/i) and a constant variance (a 2 = f3 2 / 'At). We 
know from Eq. Q that if the dot is localized onto a single 
energy state, the mean has two possible values: /io = 
or /ii = 1. Detecting a change in the mean is therefore 
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equivalent to a simple hypotheses testing problem: 



B. Bayesian solution 



H : fi 



Mo 
Mi 



(5) 



CUSUM is based on a sequential probability ratio test 
(SPRT). SPRT, in turn, is based on the concept of log- 
likelihood ratios. These are defined by 



S n = S(y?) = In 



(6) 



where p^ 0lfJLl (y) is the probability density function (pdf) 
given means liq and /ii, respectively; S is called the suf- 
ficient statistic in the parlance of statistics [13 ; note y™ 
includes all output signals (y\ . . . y n ). Further, assuming 
that the random variables Yi,...,Y n are independent and 
identically distributed (hd), 



(7) 



Sk is the sufficient statistic for random variable Uti- 
lizing the fact that is Gaussian with pdf 



we have 

Sk 



V2 



7T(J* 



Mi - Mo 



exp 



Vk 



(yk - Mo/i)^ 

2a 2 



Mi ~ Mo 

P 2 



yk 



Mo + Mi 



Mo + Mi 



At. 



(8) 



Since the bias voltage is on at the start of a pumping 
cycle, and the algorithm has to decide when to turn it 
off, CUSUM has to detect the change Mo Mi- Ignoring 
the constant pre-factor in Eq. ([8|, it is easy to see that 
S n [Eq. ([7])] has a negative drift if the dot is in its ground 
state (yk ~ Mo = 0), and a positive one when in excited 
state (yk ~ /^i = 1); a tunneling therefore should cause 
a V-shaped profile in S n . As the signal is noisy, we use 
an appropriate threshold to mitigate false positives. 

To this end, we define 



rrik — min Sn. 

l<3<k 



(9) 



We calculate rrik using a simple procedure: at each step, 
set rrik = Sk if Sk < mk-i, otherwise it retains its pre- 
vious value. In hardware, this can be done with a single 
register and a comparator. The decision rule, at each 
time step, is 



-{ 



#0 if (Sk ~ rn k ) < h, 
Hi otherwise, 



(10) 



where h is the threshold, chosen based on the parame- 
ters of the system to mitigate false positives. We stop 
pumping if dk = H\. 



The primary difference between a Bayesian and se- 
quential solution is the use of prior probabilities in the 
former case. We approximate the system described in 
Fig. [3] as a Markov chain, shown in Fig. [4] This approxi- 
mation is good because the system is assumed to operate 
in the weak coupling regime [Eq. Q], where the cav- 
ity decay rate n dominates the emitter— cavity coupling 
strength g. 



Reset 
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FIG. 4. (Color online) Markov chain approximation. We 
approximate the energy level diagram as a Markov chain. The 
individual states represent the state of the dot along with 
some number of photons in the external mode. The state of 
the dot is represented by x, and n is the number of photons 
in the external mode; x is if there is no e~ in the dot and 1 
if an e~ is present; n is a positive integer. The parameters r p 
and r e are the pumping and emission rates, respectively. We 
combine x and n into a single variable k — 2n + x; note that 
odd k states contain an e~ while even k states have no e~ in 
the dot. The system starts initially in the state k — with 
no e~ in the dot and no photon in the external mode. 

A state in the Markov chain represents the state of 
the dot and the external mode. The state of the dot is 
represented by x, and takes values or 1 depending on 
whether an electron has tunneled through to the dot or 
not, with 1 indicating the presence of an electron. The 
external mode contains n photons, where n is a non- 
negative integer. We combine the state of the dot (x) 
and cavity (n) into a single variable using, 



2n 



(11) 



Even k indicates the absence of an electron, and odd k 
the presence of an electron in the dot. The parameters 
r p and r e are the pumping and effective emission rates, 
respectively; here, r p = Q and r e = 4g 2 / k. (We ignore 
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spontaneous emission T in this approximation as its con- 
tribution is small, see Eq. ([2]).) 

Let Pk(t) be the probability of state k at time t. The 
system starts in state k = 0, so po(0) = 1 and Pk(0) = 
for all k > 0. The Markov chain in Fig.[4]can be described 
by a set of coupled differential equations: 



dpo 
dt 
dpi 

dt 
dp2 
dt 



-r p p 
r P Po ~ r e pi 
r e Pi ~ r p p 2 



(12) 



Using these initial conditions, we can solve the system of 
equations analytically; we first solve for po, then use this 
solution to solve for pi , and so on. 

Equation (13) describes the Markov chain in Fig. [4] 
when no measurement output is available. To incorpo- 
rate the information from the continuous measurements, 
we first consider a fixed time step At, and then go to 
the continuum limit. During each At the system evolves 
according to Eq. (13); then at the end of the interval 



At we update the probabilities by conditioning on the 
measured output. We now derive update formulas (con- 
ditional probabilities) for p^s using the Bayes rule. 
Let 

y(t + At) = (V x )(t + At) At + AWt+At, (13) 

where W t +At is a Wiener process with E[AWj+At] = 
and E[AW 2 +At ] = At; (P x )(t + At) is the quantum 
expectation of the state of the dot. We define 

(14) 
(15) 



?0 P(X = 0) = po + P2 

h =p(x = 1) = Pi + Ps 



where Qi and Qo represent the probabilities of an elec- 
tron to be present in the dot or not (corresponding to x 
being 1 or 0). They obey the simple equations 



(16) 
(17) 



dQo = (- r p Qo + r e Qi) dt, 
dQi = (+r p Q - T e Q x )dt, 

with initial conditions Qo(0) = 1 and Qi(0) = 0. We 
have used Eq. (13) in the above derivation. 



Since y is Gaussian random variable with means 0/At, 
we have 



p(y\x = 0) = 
p(y\x = 1) = 



1 



V27RT 2 



exp{V/2<r 2 }, 



(18) 



1 



V27RT 2 



exp{-(^-At) 2 /2a 2 }, (19) 



where, a 2 = E[y 2 ] ~ f3 2 At is its variance. Also, 

p(y) — p(y\x = 0)p(x = 0) + p(y\x = l)p(x = 1) 



1 



V2 



(Q exp{V/2a 2 } + 

Q ieX p{-(y- At) 2 /2a 2 }), (20) 



so that 

p(y\x 0) 
p(y) 



Qo + Qi exp 



y 2 - (y - At) 2 
2a 2 



which using the relation a 2 = f3 2 At becomes 



Qo + Qiexp^ Q2 



} exp {^} 

Assuming that At <c /? 2 , we get 
» Qo + Qi ( l + + ^ I I l 



i -i 



y_ 

P 2 



Qi 



At 

w 

0{M 3 ' 2 ) 



At 



0(At 3 / 2 ) 



1 - Qi(y - Q\At) /p 2 . 



(21) 



In the above derivation, we used the fact that the total 
probability is conserved, and hence Qo + Qi = 1 at all 
times. Similarly, 



p(y\x = 1) 
p(y) 



1 + Qo(y - (l-Qo)At)//3 2 . (22) 



Therefore the update formula for Qo is 

'p(y\x = {)y 

Jo 



and that of Qi is 



p(y) 

o - QoQi(y - QiAt)/f3 2 



(23) 



p(y\x = 1) 



Qi 



p(y) 

h QoQi(y - (i-Q )At)//3 2 . 



(24) 



We now combine the system dynamics given by 
Eq. (17) with the update formulae due to observation 



from Eqs. (23) and (24) to get 



Qi(t + At) = Qi(t) - (r p + r e )Qi(t) At + r p At 

+ 0.W(l-0.(«))( »" ) -£ lit)M )- (25) 

We can infer Qo straightforwardly using Qo + Qi = 1. 
Since x only takes values and 1, its expected value x is 



E[x] = lxp(x = l) =Qi. 



(26) 



We can now make precise the question of when to stop 
pumping. We need to calculate the individual probabili- 
ties pk including the prior probabilities [Eq. (|l3)] and the 
update formula due to continuous monitoring of the dot. 
Here, we use (n, x) notation instead of k to represent a 
state in the Markov chain (Fig. |4|, for clarity. Applying 
the Bayes rule we get 

p(y\n,x) p(y\x) 
p(n,x\y) = p(n,x) = —/—p(n,x). (27) 



p(y) 



p(y) 
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We have used the fact that the observation y is inde- 
pendent of number of photons n in the external mode. 
Switching back to the variable k we get 



p(k\y) 



p(y\x) 
p(y) 



p(k). 



(28) 



Using the above relation with update formulas in 
Eqs. ( |21[) a nd (22), along with prior probabilities given 
in Eq' (|13|), we get 



Po (t + At) = p (t) - r pPo (t)At 
'x{t) (y - xAt) 



■Po(t) 



pi(t + At) = pi(t) - (r e pi(t) - r p p (t))At 

'(i - x(t))(y(t) - x(t)At) 



-Pl(t) 



P 2 



p 2 (t + At) = p 2 (t) - (r p p 2 (t) - r e pi(t))At 

'x(t)(j/(t) - x(t)Aty 



-p 2 (t) 



/? 2 



(29) 



(30) 



(31) 



Here, /3 represents the quality of measurement, effectively 
the inverse of the signal-to- noise ratio (SNR): the higher 
the value of /?, the lower the signal quality. Note that as 
/3 is reduced (better signal quality), the terms from the 
(Bayesian) update formula dominates, indicate that our 
estimate of the probabilities p^ derives mostly from our 
observation; while as j3 — >• oo the probability pk converges 
to its a priori solution Eq. (13). Therefore, the Bayesian 



technique should perform no worse than the a priori solu- 
tion, at least as long as the Markov chain approximation 
remains good. 

Observe that states k G (0,1,2) can emit at most 1 
photon, while states k > 2 emit 2 or more photons. If 
we bound the tolerable multi-photon probability by e, 
then our decision circuit becomes straight-forward: we 
continue to pump until p(k > 2) = p% + p 4 + . . . < e, 
and stop pumping as soon as the inequality is violated. 
The above condition can be re- written in more convenient 
form: 



p(k < 2) = po + pi + p 2 > (1 - e), 



(32) 



and as before, we stop pumping when the inequality is 
violated. In this form it suffices to keep track of just 3 
probabilities — Po:Pi:P2 — together with the expectation 
value x. 



V. RESULTS 

In Sec. |IV[ we presented two different approaches — 
sequential and Bayesian — to improve single-photon prob- 
ability in the presence of continuous monitoring. The 
Bayesian solution (Sec. IV B) is computationally expen- 
sive, but provides a smooth transition from low-noise 



limit to the high-noise one; since in the high-noise limit 
the Bayesian updates converges to an a priori solution, 
we will do no worse than the deterministic solution 
(whe re no measurement is done). Sequential CUSUM 
(Sec. |IV A ) on the other hand, is simple and requires 
less computational resources (an integrator, a compara- 
tor and 2 registers). This is important because the deci- 
sion to turn off the pumping has to be taken in real time. 
Though the Bayesian solution is very useful for the in- 
sight it provides, particularly in high-noise settings, in 
this section we explore sequential CUSUM, as it is far 
easier to implement in real time. At high SNR it should 
approach the performance of the Bayesian solution. 

We first discuss the numerical values of the system 
parameters. Then we consider the benchmark against 
which we compare the CUSUM technique: the a pri- 
ori evolution of Eq. ( 13 ) without continuous monitoring 
(Fig.|5j. We compare the results of CUSUM to this deter- 
ministic solution, and plot the single-photon probability 
as a function of pumping rate ft (Fig. [6}; we consider 
3 cases in CUSUM corresponding to different measure- 
ment quality (/?) regimes: (i) low-noise, (ii) intermediate- 
noise, and (Hi) high-noise. The monitoring efficiency n 
has a strong effect on on the performance of the CUSUM 
algorithm, and we explore this dependence. For low n 
(= 0.1), CUSUM is not useful, and in fact is detrimen- 
tal. In this regime, the measurement output is so noisy 
that CUSUM cannot recognize when the dot becomes ex- 
cited. In such a case, increasing the monitoring strength 
(and hence the dephasing rate) nominally (and thereby 
improving 77) can lead to regimes where the technique 
performs better (Fig. [7]). Note that in all our simula- 
tions, we use the constraint 



p(2+) < 1%, 



(33) 



That is, we require that all solutions satisfy the condition 
that the multi-photon probability can at most be 1%. We 
use the fourth-order Runge-Kutta integrator rk4 [14] to 
numerically integrate the stochastic master equation ([!]). 



A. Parameter values 

In our model, five parameters characterize the system: 
g, ft, 7, T, and ft. We rescale all the parameters with 
respect to the cavity decay rate ft, which establishes the 
dimensionless (frequency) units for the simulation. In 
these dimensionless units, the parameter values are: 
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0.1, T = 0.001 and ft = 1.0. 



(34) 



The value of ft in physical (frequency) units is ~ 95 KHz 
(in [IS], / = 220 MHz and Q = 2300; in frequency units, 
ft = f/Q) . The physical units of other parameters can 
be obtained straightforwardly by rescaling with respect 
to ft. 

Pumping rate ft and measurement strength 7 are 
interesting from our standpoint, in that the effective- 
ness of CUSUM as compared to an a priori strategy 
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is strongly affected by them; Q determines the time- 
window for the decision circuit while measurement qual- 
ity (/3 = (777) -1 / 2 ) influences the ability to make the right 
decision, that is, to turn off the pumping at the right 
time. The a priori solution also improves with higher 
further eroding the benefits of the CUSUM protocol. 



B. Deterministic solution 




FIG. 5. (Color online) Deterministic evolution. We plot 
photon probabilities — p(0), p(l) and p(2+) — as a function 
of the pumping duration; the parameter values are given in 
Eq. (34) and Q = 0.1. The system is initially assumed to 
be in the state \G, 0,0). Initially p(0) = 1 and decreases as 
pumping duration is increased; p(l), on the other hand, starts 
at 0, increases initially, reaches a maximum and then starts to 
decrease; p(2+) increases with pumping duration (more grad- 
ually than p(l)) and starts to dominate at long times. The 
maximum of p(l) is « 0.73 and it does so at t ~ 19.5; how- 
ever, p(2+) is about 12%, which is quite high. Imposing the 
constraint in Eq. (33), the best time to stop pumping at these 
parameter values is t ~ 8, where p(l) « 53% and p(0) « 46%. 

To evolve the system deterministically, we integrate 
E[dp] where dp is defined in Eq. 0; due to the ex- 
pectation E[.], the stochastic contributions vanish and 
the equation reduces to the usual deterministic Lindblad 
master equation [TP] . 

Figure [5] shows the photon number probabilities — p(0), 
p(l) and p(2+) — as a function of pumping duration for 
parameters values defined in Eq. ( [34] ) and Q = 0.1. 
The probabilities are obtained by integrating the Lind- 
blad master equation with pumping turned on up to 
time t, and continuing the simulation for a sufficiently 
long time after the pumping is turned off. The initial 
state is \G, 0,0). We see that p(0) falls exponentially 
as the pumping time is increased, while p(l) increases 
initially, reaches a maximum, and then falls off; p(2+) 
increases slowly but grows to 1 at long times. The max- 



imum of p(l) is about 73% at t ~ 20; however, p(2+) 
is about 12%, which is unacceptably high. Imposing 
the constraint Eq. (33), we find that the best solution 
is p(l) ~ 53%, and is achieved for a pumping time t ~ 8. 



C. CUSUM performance 

We analyze the performance of CUSUM by integrat- 
ing the stochastic master equation ([T ) for the parameter 
values in Eq. (34), with Q £ [0.01, 0.1]; we set 77 = 1, and 
impose Eq. (33). We assume that the system is initially 



decoupled and starts in |G, 0,0). The decision rule given 
in Eq. (10) is used to stop pumping. This rule requires 



us to specify the threshold value h : to avoid false posi- 
tives; we find a close-to-optimal value for h by numerical 
exploration. For each we continue to integrate each 
trajectory for sufficiently long time after the decision cir- 
cuit stops pumping to calculate the photon probabilities. 
We repeated this procedure for 1000 trajectories and av- 
eraged to obtain pi(Q). We plot this in Figure [6] 

In Fig.[6j we consider 4 cases corresponding to different 
measurement quality /3 = (777) -1 / 2 : (i) deterministic (no 
measurement), (ii) low-noise 7 = 10, (Hi) intermediate 
noise 7 = 1, and (iv) high- noise 7 = 0.1. The better 
the quality of measurement (lower /3), the more accurate 
is our knowledge of the state of the dot; naturally, we 
expect CUSUM to perform well for lower f3 and to fare 
badly when (3 is high. This is borne out in Fig. [6] where we 
observe case (ii) performing the best and (iv) the worst. 
In fact, (iv) does worse than the deterministic case; this 
means that we are better-off not using CUSUM in the 
high- noise limit. 

Observe that cases (ii) and (Hi) provides greater im- 
provement for lower fi, and the improvement reduces as 
Q is increased. Physically, a smaller pumping rate means 
that the successive e~ tunneling events are spread out 
in time, giving sufficient time for the decision circuit to 
make the right decision; higher Q means less time be- 
tween successive e~ tunneling events, and consequently 
a tighter time- window for the decision circuit. Also, as 
we have discussed above, for high Q the optimal strategy 
is strong pumping for a short time, and we do not expect 
monitoring and feedback to yield much improvement. 



D. Measurement Efficiency 

The results in the previous section assumed that the 
monitoring efficiency was 77 = 1. This is unrealistic even 
in principle, because any realistic quantum dot will have 
multiple sources of dephasing, and we cannot have access 
to information that is lost to the dot's internal degrees 
of freedom. On top of this difficulty in principle, efficient 
monitoring is hard to do in practice. An obvious question 
is: how does lower v\ affect performance? 

To understand the effect of 77 we simulated 4 cases: 
(i) deterministic (7 = 1.0, 77 = 0), (ii) low monitoring 
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0.01 0.02 



0.04 0.05 0.06 0.07 
Q (pumping rate) 



FIG. 6. (Color online) CUSUM performance. We plot the best p(l) given constraint Eq. ( |33| for different pumping rates 
Q £ [0.01,0.1]; we find a close-to-optimal value for the threshold h by numerical exploration. There are 4 cases corresponding 
to different measurement qualities f3 — (7?7) _1//2 (here rj = 1): (i) deterministic (no measurement), (ii) low-noise (7 = 10), 
(Hi) intermediate noise (7 = 1.0), and (iv) high- noise (7 = 0.1); we plot the deterministic case for comparison. Low- noise 
(ii) naturally leads to the best performance while high-noise (iv) performs the worst; in fact, (iv) performs worse than the 
deterministic case (i); in the intermediate regime (Hi), CUSUM does better than cases (i) and (iv). In cases (ii) and (Hi), 
the performance improvement is higher for lower Q and the improvement reduces as Q is increased. Physically, lower pumping 
rate means that the successive e~ tunneling events are spread out in time, thus giving the CUSUM decision circuit has longer 
time to make the right decision; for higher rates, this time-window for CUSUM decision circuit is diminished, leading to 
lesser performance improvement. Also, the deterministic solution is better for higher Q, leaving less room for improvement by 
CUSUM. 



strength but high efficiency (7 = 1.0, 77 = 1.0), (Hi) low 
monitoring strength and low efficiency (7 = 1.0, 77 = 0.1), 
and (iv) higher monitoring strength and moderate effi- 
ciency (7 = 2.0,7^ = 0.5). Case (i) represents the de- 
terministic evolution, where we do not monitor the dot, 
though some natural dephasing occurs (Sec. VB). Cases 
(ii), (Hi) and (iv) involve continuous monitor ing a nd im- 
plement the CUSUM-based technique (Sec. VC). Case 
(ii) is the scenario of high monitoring efficiency (which, 
as we argue above, is likely unrealistic in practice), while 
in case (Hi) the efficiency is low. Finally, in case (iv) the 
monitoring efficiency is moderate, but we have boosted 
the measurement strength to compensate. 

Physically, dephasing is intrinsic in any realistic sys- 
tem, as information about the state of the dot is lost 
to various internal modes, such as phonon modes. In 
addition, there will be additional dephasing due to the 
presence of monitoring. In these terms, case (i) has in- 
trinsic dephasing but no monitoring; case (ii) has effi- 
cient monitoring and no intrinsic dephasing; case (Hi) 
has monitoring with low efficiency, but no additional in- 
trinsic dephasing. In case (iv) we assume that there is 
both monitoring and intrinsic dephasing, so the total de- 



phasing rate 7 is higher than in case (ii), and the de- 
tector efficiency 77 = 0.5 is moderate, reflecting the fact 
that only information from the monitoring is available, 
and not the information lost to the internal degrees of 
freedom. 

Figure [7] displays the single-photon probability as a 
function of Q for these four cases. The parameter values 
are defined in Eq. (34), and we have set g = 0.1. This 



plot explores just a representative slice of the (g, Q) space 
that we considered earlier. We see that case (i) performs 
consistently better than (Hi), while (ii) unsurprisingly 
outperforms both (i) and (Hi). However, case (iv) per- 
forms even better, and consistently dominates (ii). 



VI. CONCLUSIONS 

CUSUM is a simple yet powerful method that can sig- 
nificantly improve single-photon probability using contin- 
uous monitoring. The protocol itself requires only sim- 
ple components like an integrator and subtractor, along 
with 2 registers and a comparator. This technique is 
useful in the weak coupling regime when the pumping 
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rate is comparable to coupling strength. In regions of 
strong coupling, or strong pumping, it is ineffective and 
can be worse than no monitoring at all. We included var- 
ious decoherence processes in our simulations, including 
spontaneous emission, dephasing and cavity leakage. We 
modeled the electrical pumping and continuous monitor- 
ing as a stochastic master equation. Numerical simula- 
tions showed that CUSUM performs quite well in regions 
of low and significant improvements in single-photon 
probability were observed. We also studied the effect of 
imperfect monitoring efficiency on performance. 



0.02 0.03 0.04 0.05 0.06 0.07 
Q (pumping rate) 



FIG. 7. (Color online) Monitoring efficiency and per- 
formance. In this plot, we explore the effect of moni- 
toring efficiency on performance. The 4 cases considered 
are: (i) (7 = 1.0, rj = 0), (ii) (7 = 1.0, rj = 1.0), (Hi) 
( 7 = 1.0, 77 = 0.1), and (iv) (7 = 2.0, 77 = 0.5). Case (i) 
performs consistently better than (Hi) while (ii) dominates 
both (i) and (Hi). However, case (iv) performs even mod- 
estly better than (ii). Hence, a small increase in monitoring 
strength can compensate for lower efficiency. 
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